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Abstract 

We investigate the properties of the ground state of a system of in- 
teracting bosons on regular lattices with coordination number k > 2. 
The interaction is a pure, infinite, on-site repulsion. Our concern is 
to give an improved upper bound on the ground state energy per site. 
For a density p a trivial upper bound is known to be —kp{l — p). We 
obtain a smaller variational bound within a reasonably large family of 
trial functions. The estimates make use of a large deviation principle 
for the energy of the Ising model on the same lattice. 



PACS: 05.30.Jp, 05.50.+q, 75.10.Jm 



1 Introduction 



Bosons on a lattice interacting via an infinite on-site repulsion (hard-core 
bosons) represent a system of double interest. It is the simplest example 
of an interacting Bose gas and, thus, the most promising candidate for a 
rigorous treatment of Bose-Einstein condensation of interacting particles. On 
the other hand, the model is known to be equivalent to a system of ^ spins 
^ coupled via the XY- and possibly the Z-components of neighboring spins 
and exposed to an external magnetic field in the Z direction. Ordering of the 
planar component of the spins is equivalent to Bose-Einstein condensation 
or the appearance of off-diagonal long-range order (ODLRO) in the system 
of bosons. In spite of a long and extensive study the results about ordering 
are far from being complete. Apart from some exceptions, like bounds on 
the density of the condensate p[ or the discussion of the model on the full 
graph P, 1^, the most interesting and difficult results were formulated in 
spin terminology |^. These works made use of a particular symmetry, 
the reflexion positivity. This introduced severe limitations as to the value 
of the external field (zero field) and the lattice type (essentially hypercubic 
lattices). Translated into the language of the boson gas, ODLRO was shown 
at half-filling on hypercubic lattices in the ground state in and above two 
dimensions, and for low enough temperatures above two dimensions. 

In this paper we apply the boson terminology. Let L be an infinite lat- 
tice which, for the sake of simplicity, will be supposed to be regular with a 
constant coordination number (valency) k. Throughout the paper A denotes 
a finite part of L equipped with a periodic boundary condition so as to keep 
the valency constant (not really essential). The Hamiltonian we are going to 
study is 

H = - E (Kby + KK) (1) 

{x,y}eEA 

We write x,y, . . . for the vertices of L, and EA for the set of edges of A; 6* 
and create, resp., annihilate a hard-core boson at x. Boson operators at 
different sites commute with each other while 



b:b,+bj: = i (2) 

accounts for the hard-core condition. Correspondence with spin models is 
obtained by setting b^ = S~ and 6* = S^. The Hamiltonian conserves the 
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number of bosons, 

N=j:n. = Y: hih^ (3) 

and is also invariant under particle-hole transformation. We can, therefore, 
fix so that p = A^/|A| is between and |. (Here and below, if A is a finite 
set, |y4| denotes the number of its elements.) Our concern in this paper is to 
provide nontrivial upper bounds to the ground state energy, Eq. That such 
an estimate may be useful in the study of qualitative properties of the ground 
state, was an interesting point of 0. 

Let X,Y, . . . denote A^-point subsets of A, called also configurations. A 
convenient basis is formed by the states 

0(X) = n &:^vac (4) 

where \Efvac is the vacuum state. Variational estimates of the ground state 
energy are of the form 

E, < {^P\H\^)/{i;\i;) . (5) 

A trivial choice is 

V^ = ^0(X) (6) 

where the summation goes over all A^-point subsets of A. It yields (cf. Sec- 
tion 2) 

Eo < -k\A\ (^n) = ~ P^^^^ + ^^^^ • 

The bound is nothing else than minus the average size of the boundaries of 
the configurations: We call the boundary of X the set of half-filled edges and 
denote it by dX. Hence, is equivalent to 

\Eo\>WX\^i^^) (8) 

with summation over A^-point configurations. If A is a full graph, i.e., any two 
sites are neighbors, is an exact eigenvector and (H) holds with equality. 
In any other case we have a strict inequality. 

The present work is about trying to improve on this bound, i.e., to make 
the right-hand side of (§) larger by an amount proportional to |A|. It is 
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important to know that this is possible for any k and any p. In the opposite 
case, if —kp{l — p) were to be the true ground state energy per site, one could 
easily conclude that as in the full graph, the Hamiltonian (|1|) has a product 
ground state in infinite volume, 

^ = IKv^ 1^- = 1) + V^W \n. = 0)) (9) 

X 

with ODLRO and the value of the order parameter at its theoretical maxi- 
mum (p(l — p), cf. 0]). 

Since by putting equal weights on every configuration, as in (H), yields 
the average boundary size (|[), we expect that a larger value can be obtained 
by giving larger weights to configurations with larger boundaries. Therefore, 
we calculate variational bounds by using trial functions of the form 

^P, = Y.'"xci>{X) vx = v{\dX\) (10) 

i.e., vx depending on \dX\ only, and v{n) being concentrated on > \dX\. 
Computability depends on our ability to estimate the number of configura- 
tions with a given size of boundary. The logarithm of this number is the 
entropy of the Ising model in a microcanonical ensemble with a fixed mag- 
netization, Y.x£k^x = 1^1(1 ~ 2p). Indeed, if we put cr^: = — 1 for x in X 
and (Ta; = 1 elsewhere, we get an Ising configuration with the corresponding 
union of contours dX whose total length, \dX\, is the energy of the Ising 
configuration. In mathematical terms, the distribution of \dX\ satisfies a 
large deviation principle whose rate function is, apart from a shift, minus 
the specific entropy of the Ising model, cf. 0, H, The exact form of the 
entropy is unknown for two- and higher- dimensional lattices. To circumvent 
this problem, we use an approximate formula for the probability of having a 
boundary of length n, 

P^A\dX\ =n)^ exp{- ^^~jf^' } . (11) 

Here Z is for normalization and 

M = ]dX\ = kp{l - p)\A\ = {\dX\-\dX\y . (12) 

What is the approximation in the above expression? First, we replace the 
smooth and concave specific entropy s(e), having a maximum at the specific 
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Ising energy = kp{l — p), by a parabola. Since the improved bound on 
|i?o|/|^| we are going to derive is (1 + S)em with a S never exceeding 0.2, 
cf. Table 4 below, this seems to be a consistent approximation. Second, we 
surmise that the second derivative of s(e) at the maximum is —\A\/D'^. An 
analogous statement holds true for the large deviations of a sum of identically 
distributed independent random variables. Now \dX\ is the sum of identi- 
cally distributed random variables, although not independent but finitely 
dependent ones, see Section 3. So this approximation is hopefully good. In 
any case, we use the formula (|lT]) without looking for further justification. It 
may happen that energies of trial states depending on \dX\ only can be quite 
close to Eq. Then the error of the approximation ([TTD may invalidate our 
estimates as rigorous upper bounds. They still can be useful as approximate 
formulas showing the k- and p-dependence of the ground state energy. We 
apply the term 'bound' with this reservation. 

In Section 2 we present the general setup for the estimates. In Section 3 we 
derive an expression for which seems to be valid for any vertex-transitive 
graph (whose all vertices are equivalent). This formula may also be of inter- 
est in site percolation problems and in the approximation of the entropy of 
the Ising model in nonvanishing magnetic fields. In Section 4 we compute 
variational bounds by using for v{n) step-, exponential and Gaussian func- 
tions. The last two will be seen to yield improved bounds for any value of k 
and p. Section 5 summarizes the results and indicates the way they can be 
extended to the grand-canonical ensemble. 

In the true ground state vx is not a function of \dX\ only. The detailed 
dependence on X may not be relevant for computing the ground state energy 
but is absolutely crucial for qualitative properties, like ordering. Wave func- 
tions of the kind (|l^) trivially show an off-diagonal long-range order. In the 
true ground state fluctuations around a function of |i9X| destroy ODLRO in 
one dimension and may decrease the order parameter in higher dimensions. 
A discussion of the ground state wave function will be given in a subsequent 
paper 
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2 Setup for the energy estimates 



It seems advantageous to consider the eigenvalue problem of (|l]) as a problem 
of graph theory. The A^-point configurations are vertices of a huge graph 
that we call the power graph of A of order and denote by G = Ga^at. Two 
configurations X and Y form an edge of G if F can be obtained from X by 
moving a single particle from a site of X to a neighboring site, unoccupied 
in X. Thus, X and Y have A^ — 1 common sites and differ on an edge of A. 
If VG and EG denote, respectively, the set of vertices and edges of G then 

The boundary dX of X can be given a new interpretation as the set of 
neighbors of X in G. So if X and Y form an edge then Y G dX and 
X G dY . The action of —H on G is that of the usual lattice Laplacian with 
the exception that there is no subtracted diagonal term. The matrix of —H 
in the basis (^) is the adjacency matrix A = A{G) of G, that is, Axy = 1 if 
X and Y are neighbors and zero otherwise. We are interested in the largest 
eigenvalue of A, Ai = \Eo\ and in the corresponding eigenvector, a = (ax)- 
G is connected if A is connected (that we suppose), therefore A is irreducible 
(ergodic) and the Perron- Frobenius theorem applies: Ai is nondegenerate and 
largest to absolute value, and ax > for all X. We note that G is bipartite 
if and only if A is bipartite, and hence — Ai is an eigenvalue of A if and only 
if A is bipartite. 

We call nmin, resp., rimax the minimum, resp., maximum of \dX\ among 
the A^-point configurations. Clearly, n^in = 0{N^'''~^^^'^) if L is a d-dimen- 
sional lattice, and (p < |) 

nmax <kN = kp\A\ . (14) 
We have the trivial inequalities 



< Ai < n^ax . (15) 
The first is the variational bound (^: Setting f x = 1 we find 



\dX\ (16) 
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which is the same as 



{v,v) \VG\ ^^^.^^^^^^^ \VG\ ■ ^ ^ 

The upper bound in (|15D follows from the eigenvalue equation: Let Xq be a 
configuration on which ax reaches its maximum, then 

^i^Xo = ^ ay < |<9Xo|axo • (18) 

YedXo 

We note that for p = | there is a better upper bound, Ai < i(A; + 1)|A| (0, 
Theorem C.l). 

In Section 4 we shall see that nmax plays a role in optimizing the lower 
bound to Ai. It is therefore important to know rimax exactly. In there 
is equality if p is small enough. In particular, rimax = kp\A\ for all p < | 
on bipartite lattices, and it is an easy graphical exercise to see that equality 
holds for p < I on the triangular {k = 6) and on the Kagome {k = 4) 
lattices. Moreover, for both lattices ^max is constant between the densities 
^ and |: 2|A| for the triangular and ||A| for the Kagome lattice. This can 
be seen from the following argument. In general, —nmax is the ground state 
energy of the antiferromagnetic Ising model under the restriction that the 
magnetization is fixed, J^xeA^^x = (1 — 2p)|A|. However, we do not need to 
deal with the restriction. In both cases the (unrestricted) ground state is 
known to be highly degenerate. Among the exponentially large number of 
ground state configurations there are nonmagnetized ones, corresponding to 
p = |, others with concentration of down-spins p = |, and between these 
two limits p can vary by steps of -j^. The rule is to flip zero-energy spins one 
by one. The common energy of all these configurations is easy to compute 
from the fact that in each triangle there is precisely one unsatisfied bond. 
This fixes the value of nmax as given above. 

The variational bound reads Ai > B{v) where 

5(-)-%^ = E-x E (19) 

v^' ^/ X YedX X 

In the estimations below a crucial role is played by the inequality 

\\dX\-\dY\\<2{k-l) if YedX. (20) 
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When passing from X to y a neighboring particle-hole pair is interchanged. 
For both the particle and the hole the number of neighbors of the opposite 
kind can change by at most k — 1, whence (^) follows. If vx = v^dXl) and 
v{n) is a nondecreasing sequence then 



Z^n=ran,in 




f 1 








max 
— '^min 




v{ny 



where f2„ denotes the set of configurations with boundary length n. 

Above it is understood that v{n) = v{nmin) if n < n^am- We have com- 
puted bounds given by the right member of (pT]), using step-functions and 
functions with an exponential or a faster increase. In all cases we have found 
no improvement with respect to the trivial bound if p was in a neighborhood 
of |. Apparently, we have lost too much in the inequality (pT]). 

There is, however, a way to compute B{v) in leading order (|A|) by making 
a further hypothesis on some details of the large deviation principle (|ll]). 
Let z/j(X) be the number of those neighbors of X having a boundary length 
\dX\ + i. (So I'iiX) = \dX\.) For vx = v{\dX\) we have 

""^"^ = llnvinym • 

Equation (|TT]) is equivalent to 

l^nl/l^ml ~ exp{-[(n - Mf - (m - Mf]/2D'^} . (23) 

An analogous statement for the ratio of z/j(X) and i^jiX) is certainly wrong 
for all X separately but may be correct for the ratio of their sums over Qn- 
In this hope we formulate the hypothesis that 



2k + 2<i,j <2k-2 . (24) 



e 20-= 

-1 



The third member of ( ^41) is obtained by dropping a term of order |A 
With (0) or (^) and (^ Eq. (|22|) reads 

B{v) = 2D^v{n) /^'^e 202 v{n){v)n (25) 




n-M • 



= ^ e^'v{n-i). (26) 

This is the starting point of the estimates of Section 4. 
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3 Statistics of the boundary lengths 



In this section we show that the mean square deviation of \dX\ is given by 

D"^ = [k- 2{2k - l)p(l - p)]kp{l - p)\A\ = [k- 2{2k - l)p(l - p)]M . (27) 

We have found this expression first for ci- dimensional hypercubic lattices, 
and checked it later to hold also for the triangular, honeycomb and Kagome 
lattices. This is somewhat surprising because our derivation below needs 
knowledge of the rather different local neighborhoods up to next-nearest 
neighbors. The only common feature of all these lattices seems to be that 
all sites are symmetry-related and thus equivalent. Therefore, it should be 
possible to prove Eq. (^) on the basis of vertex-transitivity alone. 

Equation ( p7|) is obtained by using the grand-canonical probabilities Pk,p^ 
i.e., by filling the sites of A independently and with equal probability p. We 
expect smaller order corrections to appear if the canonical distribution is 
used. So in this section X is a random subset of A whose probability to be 
selected is — p)'"^'"'"^', and rix = nx{X) is a random variable taking the 

value 1 if X is in X and otherwise. Then all are independent and take 
1 with probability p and with probability 1 — p. We define 

/x = ^ (1 - uy) (28) 

y&dx 

where dx denotes the set of neighbors of x in A. Clearly, 

U = kp{l - p) . (29) 
The boundary length of X is obtained as 

\dx\ = y: MX) . (30) 

Thus the mean value of f|30f) is 



M = = kp{l - p)\A\ (31) 

xeA 

as found earlier. 

Let d{x, y) denote the graph distance of x and y in A, i.e., the length of the 
shortest walk between them. Since fx and fy are independent if d{x, y) > 2, 
we find 

D'=j:if^-Tx)ify-Ty)= E (32) 

x,y:d{x,y)<2 
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r{x, y) = fxfy - fx ■ (33) 

The computation of the different terms is straightforward by observing 
that nl = Tlx, (1 — nxY = 1 — Ux and nx{l — Ux) = 0. The contribution of 
the diagonal terms x — y is the same for any /c-regular lattice. Namely, 

fl=kp{l^p) + k{k-l)p{l-pf (34) 

J2 r{x, x) = \A\r{x, x) ^ [k - {2k - l)p + kp^]M . (35) 

The contribution of nearest neighbor pairs depends on the number of triangles 
containing a given edge. If there are £ such triangles then 



= p'm - p) + P -ly-m- pY] m 

^ r{x, y) = k\K\T{x, y) = p[lp - {2k - 1)(1 - p)]M . (37) 

x,y: d{x,y)=l 

If X and y are next-nearest neighbors to each other, they may have m common 
nearest neighbors. Then 

fj-y = p^[m{l-p) + {e-m){l-pf] 
r{x,y) = mp^(l-p) . (38) 

In (i-dimensional hypercubic lattices {k = 2d) there are next-nearest neighbor 
pairs with m = 1 and m = 2. Their contribution to is 

Yj f{^^y) = ^M'r{x,y)m=\ + '^\V\\Mr{x,y)r 

x,y:d{x,y)=2 \ / 



)m=2 

= p^M + {k- 2)p'^M ^{k- l)p^M . (39) 
For the triangular lattice {k = 6) 

r{x, y) = k\K\ [r (x, y)m=i + r(x, y)m=2] ^ {k - 3)p^M . (40) 

x,y:d{x,y)=2 

In the honeycomb lattice {k — 3) each site has 6 next-nearest neighbors, all 
of the type m — 1. So 

Y r{x, y) = 2A;| A| r{x, y)m=i ^ {k - l)p'M . (41) 

x,y: d{x,y)=2 
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In the Kagome lattice {k = 4) there are 8 next-nearest neighbors with m = 1: 
^ r{x, y) = 2k\A\ r{x, y)^=i = (k - 2)p^M . (42) 

x,y:d{x,y)=2 

Finally, we obtain by adding (|35|), (0) with £ = and (|3|) for 
hypercubic lattices, (^), (^) with i = 2 and (|40|) for the triangular lattice, 
(II), (13) with £ = and (|1]) for the honeycomb lattice and (gSj), (0) with 
£ = 1 and (H) for the Kagome lattice. All yield {^). 



4 Energy estimates 

By inspecting Eq. (p5|) it is clear that v{n) has to be chosen in such a way 
that it shifts the expectation value of the Gaussian upwards. The appropriate 
choice can be either a rapidly - at least exponentially - increasing function 
or a function concentrated on values well above M. If in the numerator we 
had v{n) also at the place of {v)n, B{y) = rimax could be reached. However, 
typically {v)n/v{n) < 1 and decreases as we modify v by putting increasing 
weights to larger boundary lengths. This limits the maximum of B[v). We 
repeat here the DLS-bound 0, mentioned in Section 2, 



and the trivial upper bound (|T5|), 



(43) 



Ai < M 



l-p 



(44) 



The right members of these inequalities are upper bounds to B{v) as well. 
Step functions 

Let v{n) = 1 ii Til < n < n2 and otherwise. Suppose that rii — M = 
m > ci\A\ and n2 — > C2\f\A\ where Ci and C2 are positive constants. We 
have 



"2 



B{v) = e ^ 



v,ri=?il 



n=ni 



2D2 



n-M ■ 

71- 



n- 



(45) 



Because of the lower cutoff, the Gaussian is sharply concentrated on rii. 
Except for v{n — i), we can therefore replace n by rii in the square bracket. 
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Also, with the above choice of n2, the dependence on n2 is neghgible. We 
find after some manipulations that apart from smaller order corrections 



B{v) = Hi 



=-2fc+2 / i=l 



(n-M)2 
2D^ 



(46) 



We simplify the fraction with q-'^^I'^^^ and evaluate it to find 1 — g^*"*/-^^ 
plus a correction which vanishes as |A| oo. Inserting this into Eq. (^Bf ) we 
obtain 



+ m) 



i=-2fc+2 

^2 (A;- 1)(2A;- 1) 
'M 4 A; - 3 



m 



(47) 



From Eq. (|27| ) we see that D'^/M tends to as p goes to zero. Hence, the 
coefficient of m/D"^ is positive if p is small enough and we have the best 
(improved) bound for rii close to nmax- Then m/D"^ p/k and we are not 
in confiict with (^^. On the other hand, as p goes to |, D"^ /M tends to \ 
and the coefficient of m/D"^ in (^) becomes negative. Thus, for any k >2 
there is a neighborhood of p = | in which the maximum of B{v) is reached 
for m = 0, that is, we get no improved bound. 



Exponential trial functions 

Let v{n) = e^"" with an a; > 0. Suppose first x < 



where 



M 



2D^ 



(48) 



Plugging V into Eq. 



due to the factor e the expectation value of the 



Gaussian is shifted from M to M + 2xD < rimax both in the numerator and 
in the denominator. The new expectation value can replace n in its other 
occurrences in the numerator. Apart from smaller order corrections we get 



B(v) = 



FJx) 



= M (l + 2xD^/m) Fu{x) = MGk,p{x) (x < x^^x) , (49) 



E 



2fc-2 



-2fc+2 ' 



Y^2fc-2 2xi 



cosh(2fc — |)a; — cosh(2/c — |)a; 
cosh(4A; — |)a; — cosh(4/c — |)a; 



(50) 
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Expanding Fk{x) around we find 



Gk,pix) = (l + 2xD'^/M) [1 - (A; - l){2k - l)x^ + 0(x^)] > 1 (51) 

if X is small enough, so we have an improved bound for any k and p. 

If X > a;niax, the maximum of the shifted Gaussian is outside the range 
of summation. Therefore the distribution is sharply concentrated on a small 
neighborhood of nmax, and the computation of B{y) can be done in analogy 
to the case of the step function. The result is 



B{v) 



2k-2 



-2k+2 ' 



2k-2 
. 1=0 



2k-2 



iX > Xr, 



(52) 



Notice that B{v) varies continuously with x and both (|^) and (|52D yield 



B{v) = n^^^Fk (XmaxJ 



if X = Xn 



(53) 



In order to obtain the best bound, B{v) has to be maximized with respect 
to X. Gk,p{x) has a unique maximum at some positive x. This can be the 
location of the maximum of B{v) only when it is smaller than Xmax, and then 
the optimal B{v) is computed with it from (^U]). Table 1 shows that this is 
the case for all A; if p is not too close to 0. However, for small enough densities 
the maximum of Gk,p is attained at an x above Xmax- For these densities the 
first expression (|^) increases with x up to Xmax while the second expression 
(|5^ ) decreases with x for all densities. So the highest bound is provided by 

The conclusion is that for all densities we can obtain the best bound by 
maximizing ( ^Dj ) with respect to x under the condition that x < Xmax- 

Numerical results on the best estimates are summarized in the tables 
below. They apply to bipartite lattices whenever is given by Eq. (^7]). 
According to Table 1, they are valid also to the Kagome and triangular lat- 
tices, because the largest density for which the maximizing x is Xmax remains 
below I for every k. 
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p 


k 


2 


3 


1 


6 


8 


0.1 


0.038* 


0.026* 


0.020* 


0.014* 


0.010* 


0.2 


0.120* 


0.089* 


0.071* 


0.040 


0.027 


0.3 


0.206 


0.082 


0.047 


0.024 


0.016 


0.4 


0.167 


0.057 


0.030 


0.013 


0.008 


0.5 


0.152 


0.048 


0.023 


0.009 


0.005 



Table 1. Rates of increase, x, of the exponential trial functions, maximizing 
Gk,p- Numbers with an asterix correspond to Xmax- 



p 


k 


2 


3 


4 


6 


8 


0.1 


1.106 


1.103 


1.102 


1.099 


1.098 


0.2 


1.199 


1.157 


1.128 


1.100 


1.088 


0.3 


1.155 


1.075 


1.051 


1.033 


1.027 


0.4 


1.094 


1.035 


1.019 


1.009 


1.006 


0.5 


1.077 


1.024 


1.012 


1.004 


1.002 



Table 2. Maximal values of Gk,p{x), computed with the entries of Table 1. 



p 


k 


2 


3 


4 


6 


8 


0.1 


0.199 


0.298 


0.397 


0.594 


0.791 


0.2 


0.383 


0.555 


0.722 


1.056 


1.393 


0.3 


0.485 


0.677 


0.882 


1.302 


1.725 


0.4 


0.525 


0.745 


0.978 


1.453 


1.932 


0.5 


0.538 


0.768 


1.012 


1.507 


2.005 



Table 3. Best estimates of |£'o|/|A| obtained by multiplying kp{l — p) with 
the corresponding entry of Table 2. 
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k 


2 


3 


4 


6 


8 


p 


0.23 


0.19 


0.17 


0.16 


0.15 


X 


0.158 


0.080 


0.051 


0.031 


0.021 




1.206 


1.159 


1.143 


1.129 


1.123 



Table 4- The maximum of Gk,p{^) ^ function of p and x G [0,Xmax], 
together with the maximizing p and x. In each case the latter is Xmax- 



Gaussian trial functions 

Let f (n) = e~*^"~^^i) /^^i where Mi > M and oo with increasing 

|A| but otherwise are free parameters. The product of the two Gaussians 
gives rise to a new Gaussian with mean value and variance 

MDl + MiD' DID' 
= Dl + D' = dTTW' ■ ^ ^ 



A straightforward computation yields the same B{v), Eqs. (|4^), (|52D and 
(|53| ) as for the exponential trial function with 

1 Ml - M 

So Gaussian and exponential trial functions provide the same bound, at least 
in leading order. Moreover, for any fixed x > there is a one-parameter 
family of Gaussians giving the same result. An interesting choice is Df = 
o(|A|), e.g. Df = constyjA|. In this case M2/M1 = 1 and Dj/Dj = 1 
asymptotically. 



5 Summary and concluding remarks 

We have presented variational estimates of the ground state energy of a 
gas of hard-core bosons on regular lattices. The wave functions we have 
used depended only on the size of the boundary of the iV-point configura- 
tions. Therefore, the estimates could be based on a large deviation principle 
governing the distribution of the boundary sizes. The corresponding rate 
function is related to the entropy of the (ferromagnetic) Ising model on the 
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same lattice. We have derived a formula for the mean-square deviation of 
the boundary sizes and applied it in a quadratic approximation of the rate 
function. 

The best estimates we have found are of the form (^) where for given 
k and p the parameter x is to be determined numerically so as to maximize 
B{v). This has been illustrated in Tables 1-3. The maximum is realized 
either by an exponential function or by a one-parameter family of Gaussian 
functions. 

Extension to the grand-canonical ensemble is straightforward. Adding 
—pJ2xeA^x to the Hamiltonian resumes in adding pp\A\ to B{v). After 
maximizing with respect to x we have to maximize with respect to p. More 
precisely, if b{p) is the maximum of B{v)/\A\ with respect to x, first we extend 
b{p) symmetrically to p > | and then determine b{p) = maxp[6(p) -|- pp] 
to get an estimate of the ground state energy per site in the full bosonic 
Fock space. It is more b{p) than b{p) which is useful for the equivalent spin 
model whose Hamiltonian now contains an external magnetic field in the 
Z-direction. The function p{p) that we can find in determining b{p) is only 
an approximation of the true relationship which exists between the chemical 
potential and the density in the ground state. Table 3 suggests that p(0) = |. 
This rigorously holds true for any positive temperature; even we have a much 



stronger 'uniform density theorem', as in the Hubbard model |TT|]: Because 
of the particle-hole symmetry, 

Trra^e-^^'^-o _ Tr (1 - n^)e-^-^^=o _ 1 

is valid in any A with free boundary condition. 

In our estimates the lattice structure appears only through k, the coordi- 
nation number. Therefore we obtain the same result for the Kagome lattice 
as for the square lattice, and for the triangular lattice as for the simple cubic 
lattice. Although the energy corresponds to a nearest neighbor correlation, 
the details of the lattice geometry certainly influence the p-dependence of the 
exact ground state energy per site. The mark of the lattice could be recovered 
by using the exact Ising entropies instead of their quadratic approximants. 
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